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ABSTRACT 

Recent theoretical and experimental evidence indicates that many of the materials that are 
thought to exist in the mantles of terrestrial exoplanets will metallize and become good conductors 
of electricity at mantle pressures. This allows for strong electromagnetic coupling of the core and 
the mantle in these planets. We use a numerical dynamo model to study the effect of a metallized 
lower mantle on the dynamos of terrestrial exoplanets using several inner core sizes and mantle 
conductivities. 

We find that the addition of an electrically conducting mantle results in stronger core-mantle 
boundary fields because of the increase in magnetic field stretching. We also find that a metallized 
mantle destabilizes the dynamo resulting in less dipolar, less axisymmetric poloidal magnetic fields 
at the core-mantle boundary. The conducting mantle efficiently screens these fields to produce 
weaker surface fields. We conclude that a conducting mantle will make the detection of extrasolar 
terrestrial magnetic fields more difficult, while making the magnetic ffelds in the dynamo region 
stronger. 



Subject headings: Keywords 



Introduction 



There is evidence on every planet in our so- 
lar system, except Venus, of a dynamo-generated 
planetary magnetic ffeld either today, or at some 
point in the past. The existence and morphology 
of these dynamo-generated ffelds are strongly con- 
strained by the properties of the planetary deep 
interior. With the vast array of extrasolar rocky 
planets now being discovered, an obvious ques- 
tion is what to expect for their dynamo-generated 
magnetic ffelds. Since these exoplanets provide 
a larger range of planetary properties than those 
found in our solar system, the potential exists for 
interior dynamics not seen in our solar system. 

Detectability is an important factor in the 



study of extrasolar planetary magnetic ffelds. 
There are two ways that the magnetic ffeld of an 
extrasolar planet could be detected from Earth. 
The ffrst occurs when electrons from the stellar 
wind interact with the dynamo-generated mag- 
netic ffeld from the planet, emitting cyclotron ra- 
diation dFarreh et al.|1999||Griefimeier et al.|2QQ7 



Lecacheux 1991). The radiated power associated 



with this is 



rad 



(X5' 



0.58^ 



-1.17 



(1) 



where B is the magnetic ffeld and a is the planet- 
star distance. The constant of proportionality is 



related to the strength of the solar wind (Farrell 
eraL][T999l ). 

The second method by which extrasolar plane- 



tary magnetic fields could be detected is through 
magnetospheric interactions between a close-in 
planet and its host star. This mechanism was pro- 
posed to explain the observations of "hot spots" 
in the chromospheres of some stars, which ro- 
tate with the planetary orbit. If a planet with a 
magnetic field orbits close enough to a star, it is 
possible that the magnetic fields lines may join the 
two bodies and trap plasma in the closed field lines 



between them (Cohen et al. 2009). The presence 



of the planetary magnetic field can be detected 
indirectly through interaction of this plasma with 
the host star. 

Both of these signatures of extrasolar plane- 
tary magnetic fields are more easily observable 
when the magnetic field at the planetary surface 
is stronger. 

Many of the large number of extrasolar terres- 
trial planets which have been discovered have a 
mass greater than that of Earth ( Schneider|2012 ). 
This implies that the pressures and temperatures 
inside "super-Earths" can be considerably higher 
than in any of the terrestrial planets in our so- 
lar system ( Valencia et al.|[2006 ). This raises the 
possibility of novel material properties inside the 
deep interiors of these planets that could poten- 
tially have an important effect on a planet's geo- 
dynamics. 

One interesting possibility is the pressure- 
induced metallization of mantle materials. In the 
Earth, the mantle is largely electrically insulating, 
mainly because its main constituent (perovskite) 
is not expected to metallize until pressures far be- 
yond any which are expected to be found in the 
deep interiors of even the largest rocky exoplan- 



ets (Tsuchiya & Tsuchiya 2011). In exoplanets. 



the possibility of compositions which differ sig- 
nificantly from Earth could lead to electrically 
conducting lower mantles in even small rocky ex- 
oplanets (Ohta et al. 2012). Recent studies have 



shown that there are several common minerals 
which should metallize. These include CaSiOs 



(Tsuchiya fc Tsuchiya 2011), FeO (Ohta et al. 



2012), and AI2O3 (Nellis 2010) indicating that 



lower mantle metallization may be a phenomenon 
that is common in rocky exoplanets. 

When mantle materials are metallized a num- 
ber of their physical properties change signifi- 
cantly. A metallized mantle should have a high 
thermal conductivity as well as a high electri- 



cal conductivity because of the Wiedemann- Franz 
law, which states that the thermal conductivity 
is proportional to the electrical conductivity and 
temperature. Another concern is that the slope of 
the liquidus of the material could change, poten- 
tially leading to a liquid lower mantle. This can 
be definitively ruled out in the case of FeO, where 
experimental data show no change in the curva- 



ture of the liquidus when FeO metallizes ( Boehler 



1992). There are no experimental data on the 



melting temperatures of CaSiOs and AI2O3; how- 
ever, Schreinmaker's rule ( Zen||1966 ) implies that 
a liquidus with dT/dP < is not possible. Any 
melting requires that metallic CaSiOs and AI2O3 
behave differently than FeO, and that a fortuitous 
combination of adiabat slope and liquidus slope 
exists. We feel this situation unlikely so we hence- 
forth assume that the mantle is completely solid 
in this study. 



Chan et al. (2008) published numerical dynamo 



simulations with a conducting mantle layer that 
had conductivity which varied sinusoidally in lat- 
itude and longitude. They found that the addi- 
tion of this electrically conducting mantle layer 
could cause a previously steady dynamo to vacil- 
late and, if the conducting mantle layer was thick 
enough, could stop a dynamo from operating alto- 
gether. However, the parameters this study chose 
were for the benchmark dynamo, an intentionally 
placid dynamo solution run at unrealistic parame- 
ters that is normally used to ensure that a numeri- 
cal dynamo model is working correctly. This raises 
concerns about its relevance to planetary regimes 
since the parameters in the benchmark dynamo 
are much less realistic than those used in most 
planetary studies. 

The purposes of these studies were not to model 
the metallization of the mantle, and so they used 
lower mantle conductivities than one might ex- 
pect from a metallic mantle. Furthermore, they 
concentrated on the effects of heterogeneity in the 
mantle conductivity. Here we use a numerical dy- 
namo model running at more realistic parameters 
to study the effect of lower mantle metallization on 
the dynamos of possible extrasolar terrestrial plan- 
ets with specific interest in the observable proper- 
ties of these dynamos. 

The study of terrestrial exoplanets, especially 
terrestrial exoplanetary interiors remains under- 
constrained. For a dynamo to exist on these plan- 



ets, among the most important factors is the state 
of the mantle. The power to drive the dynamo 
of an extrasolar terrestrial planet is controlled by 
the mantle, so the efficient transfer of heat out of 
the planet is of great importance. The presence 
of plate, tectonics and vigorous mantle convection 
are both efficient ways that planets can drive dy- 
namos in their cores. Currently, mantle convection 
on these bodies is poorly understood, especially 
for large exoplanets. The field remains sharply 



divided ([Lenardic fc Crowley 2012 O^Neih fc 



Lenardic||2QQ7[ [Stein et al.||2QlH |2Q13| I Valencia 



& O'Conneh 2009^ van Heck, Tackley pTT) as 



to the likelihood of plate tectonics, and the vis- 
cosity structure ( Karato||2Qll ) of large terrestrial 
exoplanets. 

The effect of a metallized mantle layer on man- 
tle convection has been studied by I van den Berg| 
(|2010|). They found that the addition of 



et al. 



an electrically conducting layer at the the bot- 
tom of the mantle caused the bottom boundary 
layer to heat, become buoyant and rise to the up- 
per mantle. This leads to an increased heat flux 
at the core-mantle boundary (CMB) which could 
increase the power available to the dynamo, but 
shorten its lifetime. This would be a secondary ef- 
fect as the heat flux increase is less than an order of 
magnitude. There are other, less well constrained 
properties of extrasolar terrestrial planetary man- 
tles which will have a greater impact on the heat 
flux from the core than mantle metallization (e.g. 
radiogenic heating, or the presence of plate tec- 
tonics). 

2. Expected Effects of Mantle Metalliza- 
tion on the Dynamo 

An electrically conducting mantle should af- 
fect the dynamo in two ways. First, any quickly 
varying components of the magnetic field should 
be screened out by the skin effect before they 
reach the surface, weakening any observed field. 
Inside the solid mantle layer, the magnetic field 
obeys a diffusion equation, with a diffusivity equal 
to 7^ = l/(crMMo) where (Jm is the conductivity 
of the layer and jIq is the magnetic permeabil- 
ity of free space. Neglecting the spherical ge- 
ometry of the core, the magnetic field is atten- 
uated in a solid conducting mantle proportional 

to e~^V^/*^^^\ where uo is the frequency of the 



magnetic field variations at the top of the core 
and d is the thickness of the conducting man- 
tle. This implies that the thicker the mantle layer, 
the weaker the observed field. Also, higher mul- 
tipoles will be preferentially damped due to their 
more rapid time variation compared to lower mul- 
tipoles (Christensen & Tilgner 2004), while the 



dipole component of the magnetic field may not 
be greatly affected. 

This screening effect also applies to non- 
axisymmetric components of the field that are 
being advected by the background flow at the top 
of the core. For example, in dynamo simulations, 
equatorial flux spots are a common occurrence, 
and typically drift westward (Finlay fc Jackson] 



2003). From the mantle reference frame these 



spots are viewed as a time varying magnetic field 
and hence, they will be screened. In cases where 
the dynamo generated field is exceptionally steady 
and axisymmetric, the screening effect may be 
unimportant as the timescale of field change will 
be very large. 

The second feature we expect when a conduct- 
ing mantle is added to a dynamo is magnetic shear 
at the core-mantle boundary due to flux freezing 
in the mantle. The magnetic field should anchor 
itself in both the solid mantle and the convecting 
liquid outer core (Moffatt 1978|. Shear should 



then be created between the solid mantle and the 
strong zonal flows which are present in planetary 
cores. Simultaneously, the fluid in the outer core 
should feel a Lorentz force from the stretching of 
magnetic field lines anchored in the mantle. While 
the screening effect is a kinematic process that 
simply acts on time varying magnetic fields gener- 
ated by other means, the magnetic coupling effect 
requires a fully coupled modeling approach. 

The thickness of the metallized part of the man- 
tle should have a significant effect on the char- 
acter of the observable field. We ffist note that 
the screening effect and the Lorenz feedback ef- 
fect scale differently with metallized mantle thick- 
ness. Electromagnetic screening becomes more 
important as the thickness of the metallic man- 
tle layer increases, since the field is attenuated by 

distance proportional to e~^^^'^'^^\ Conversely, 
the Lorenz force acting on the fluid by the mantle 
is nearly independent of mantle thickness. This 
is evident when examining the equation for the 
torque (F) on the core by the mantle (adapting 



Gubbins fc Roberts| ( p87t ) 

r=— ff Br{rxB)r^sm0d9d(j) (2) 

f^O J J 

CMB 

The surface integral is over the CMB and does not 
involve mantle thickness. 

The thickness of the metallized mantle layer in 
a given exoplanet depends strongly on the size of 
the planet, the size of the core, and the material 
which is metallized. For example, FeO should met- 
allize at approximately 55 GPa along the Earth's 
geotherm ( Ohta et al.|2Q12 ) meaning that a planet 
significantly smaller than Earth could have an 
electrically conducting mantle layer as long as it 
had a large amount of FeO in the mantle. Con- 
versely, CaSiOs is not expected to metallize un- 
til pressures of 600 GPa, meaning that a con- 
ducting mantle layer should form in CaSiOs rich 



planets with masses greater than 5M^ (Tsuchiya 
& Tsuchiya 2011). By varying the composition, 
planet size, and core size, a conducting mantle 
layer of nearly any thickness can be achieved. Be- 
cause of this we avoid fixing a planetary radius or 
mass, as the arguments here apply equally well to 
a range of planet sizes. As the surface magnetic 
field strength is strongly dependent on planetary 
size and core-mass fraction, we discuss only the 
properties of the field at the top of the electrically 
conducting region. 

3. Numerical Model 

To study the effect of an electrically conduct- 
ing mantle on dynamo generation in planets we 
use the Kuang-Bloxham numerical dynamo model 
(Kuang & Bloxham 1999). This model has been 



successfully applied to study the magnetic field of 
many of the planets in our solar system. It nu- 
merically solves the three-dimensional, nonlinear, 
Boussinesq, magnetohydrodynamic equations in a 
convecting, rotating, spherical shell, representing 
the fluid outer core. The inner core of our model 
is solid, and conducts electricity with the same 
conductivity as the outer core. 

In the model, all the equations are solved in 
their non-dimensional forms. We use the radius 
of the core (r^) as the length scale, the magnetic 
diffusion time r = rl/r] as the time scale and the 
magnetostrophic balance intensity B = y^2ftp/ac 
(where Q is the rotation rate of the planet, p is 



the density of the core, and ac is the electrical 
conductivity of the outer core) as the magnetic 
field scale. We also use /ib^o as the temperature 
scale, where Hb is the buoyancy flux at the inner 
core boundary. The non-dimensional equations of 
the model are: 



En 



d_ 

m 



u • V ) u + zxu 
dT 



dt 



-Vp + JxB + Raer + E\/^ 

u-VT = q^S/^T (4) 

Vx(uxB) + V^B (5) 



dB 



Here B, u, and J are the magnetic field, velocity 
and current density respectively, p is the modified 
pressure, and z is the rotation axis of the system. 
Also, T is the temperature, while is the temper- 
ature perturbation. 

The nondimensional numbers in Equations ([3|- 
([5]) are the Rayleigh number (Ra), the Ekman 
number (E)^ the magnetic Ekman number {Erj), 
and the Roberts number (^h:), and are given by 



(6) 
(7) 
(8) 
(9) 



Ra = 


2nr] 


E = 


V 


2nrl 


E, - 


2nrl 




K 


Qk = 


- — 



V 



where ar is the thermal expansion coefficient, go 
is the gravitational acceleration at the CMB, u 
is the kinematic viscosity, and k, is the thermal 
diffusivity. In all our models we set Ra = 12000, 
E = 2.1125 X 10-^ Er^ = 4.225 x 10"^ and qk = 5. 

Our model differs from most numerical dynamo 
models in its treatment of the outer boundary. 
We include a solid mantle layer, the electrical 
conductivity of which can be specified arbitrarily. 
We specify the relative conductivity of the mantle 
layer with (Jmc = o-m/ctc where (Jm is the con- 
ductivity of the conducting mantle layer. 

Our models use spherical harmonics as basis 
functions in the azimuthal directions, and a com- 
bination of Chebyshev expansions and compact fi- 
nite differences in the radial direction. All mod- 
els presented here have Lmax = 56, mmax = 53 
and the number of grid points in the inner core. 



outer core, and mantle are Ni = 50, No = 64, and 
Nm = 50 respectively. For numerical reasons we 
make use of scale-dependent viscosities and diffu- 
sivities, we have applied them lightly in this study, 
using them only for L > 20 in order to minimize 
their dynamical effects. 

We model a metallized mantle with a spherical 
shell of uniform conductivity on the outside of the 
dynamo region. In all models the spherical shell 
extends from Tq to l.OTro. We also make the shell 
highly conducting, varying aMC from 1/8 to 1. 
As a control, we run a model with a relatively 
insulating mantle {(Tmc = 1/400). A schematic 
diagram of the model is shown in Figure [l] 

The strong rotational influence on convection 
in planetary cores implies that the solid inner core 



has a significant effect on the dynamo ( Heimpel 
et al.|[2QQ5{ Stanley et al.|[2QQ7 ). Since exoplanets 
can be found in different stages of their thermal 
evolution (depending on their age), we use three 
different inner core sizes in this study. As a planet 
ages, the heat lost from the core will cause the 
solid inner core to grow, so the different cases can 
be considered proxies for different stages of the life 
of the dynamo. Here we use rio = 0.35, 0.5, and 
0.70. 

4. Results 

In all cases we find that the addition of an elec- 
trically conducting mantle significantly enhances 
magnetic field generation in our models. The most 
striking example of this is in the axisymmetric az- 
imuthal {(f)) component of the field near the CMB 
(Figure |2|. Here magnetic fields anchored in both 
the solid, stationary mantle, and the fluid outer 
core are sheared by the strong zonal flows at the 
top of the core. This provides a source of magnetic 
field stretching, which strengthens the magnetic 
field in the azimuthal direction. 

We find that the strength of the poloidal com- 
ponent of the field at the CMB is increased as well 
(Figure [3|. The poloidal component of the field is 
of special interest as it is the component which 
is observable from outside the conducting region 
(the toroidal field requires poloidal currents, which 
are only present in the conducting region). In our 
models, the energy in the non-axisymmetric com- 
ponent of the poloidal field increases as the con- 
ductivity of the mantle layer increases. In all cases 



the energy in the axisymmetric component of the 
poloidal field remains approximately constant. 

This preferential increase in the non-axisymmetric 
component of the field can be explained by noting 
that the toroidal field at the CMB is predomi- 
nantly axisymmetric, due to the large zonal ax- 
isymmetric flows shearing magnetic field there. 
Because a dynamo cannot create axisymmetric 
poloidal energy by any flows acting on axisym- 



metric toroidal magnetic fields (Bullard & Gell- 
man|1954), we should expect that the axisymmet- 



ric poloidal field should not increase due to an 
increase in the axisymmetric toroidal field caused 
by a conducting mantle layer. It is possible to 
create non-axisymmetric poloidal fields via non- 
axisymmetric velocities acting on axisymmetric 
toroidal magnetic fields, so an increase in the 
strength of the axisymmetric toroidal magnetic 
field should imply an increase in the strength of 
the non-axisymmetric poloidal magnetic field. 

The axisymmetry of a dynamo is an important 
factor in the potential observability of extraso- 
lar terrestrial planetary magnetic fields. For the 
magnetic field to reach the exterior of the planet 
where it can be observed, it must first be screened 
through a conducting lower mantle. As discussed 
earlier, non-axisymmetric fields are screened more 
efficiently than axisymmetric fields. This becomes 
apparent if we plot the poloidal magnetic ener- 
gies above the conducting layer (Figure [4]). We 
see that depending on the shell thickness, the ob- 
servable field is either only marginally stronger or 
much weaker than the field in a model without a 
conducting mantle layer. 

As Tio increases there is a marked decrease in 
poloidal field strength at the top of the conduct- 
ing layer relative to the insulating case (Figure [4|. 
When considered with Figure [3J the reason for this 
becomes clear. As the liquid outer core becomes 
thinner the dynamo becomes less axisymmetric at 
the top of the dynamo region and the characteris- 
tic timescale of variation becomes shorter (Aubert] 



et al. 2009). Both of these effects cause the field 



to be weakened by the screening effect of the con- 
ducting layer. This means that thinner shells are 
more susceptible to the screening effect discussed 
earlier. The increase in non- axisymmetry with in- 
creasing rio has been observed previously in stud- 
ies which modelled the evolution of the Earth's dy- 
namo through time (jAubert et al. 2009 Roberts 



fc Glatzmaie7][2QQT] ). 
5. Conclusions 



We have used a numerical planetary dynamo 
model to investigate the effect of the metalliza- 
tion of silicate mantles on the observable magnetic 
fields of super-Earths. We have carried out mod- 
els for three different inner core sizes to simulate 
these planets at different stages in their thermal 
evolution. In all cases we find that the strength of 
the internal magnetic field increases substantially, 
owing to the magnetic shear provided at the CMB 
by the conducting mantle. We also find that the 
addition of a conducting mantle makes the field 
significantly less axisymmetric at the top of the 
dynamo region. After being screened through the 
conducting mantle layer we find that the observ- 
able field shows either a modest increase in field 
strength (at r^o = 0.35) or a significant decrease 
in field strength (at Vio = 0.7). 

As we have used a thin conducting layer in 
these models (compared to the range that is 
possible for terrestrial exoplanets), we expect 
that in larger planets, the screening effect would 
be even stronger than we observe here. This 
means that any planets with a metallized man- 
tle should have surface fields which have been 
significantly weakened by a combination of the 
non-axisymmetrization of the dynamo, and the 
screening effect of the mantle. We therefore ex- 
pect that the metallization of silicates should make 
the detection of dynamo-generated magnetic fields 
from super-Earth's more difficult than previously 
anticipated ( Driscoll fc Olson||2011 ). 
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Fig. 1. — Schematic diagram of the structure of 
a planet with a conducting mantle layer. Our nu- 
merical model solves in the region below tm- 
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Fig. 2. — Contours of the axisymmetric azimuthal 
{(p) component of the magnetic field in the core for 
an insulating mantle model {(Tmc = 1/400, left) 
and for a conducting mantle model {ctmc = 1, 
right) at a single moment in time. The shaded 
regions indicate the inner core (center) and the 
conducting mantle layer (outside). Note the dif- 
ference in scales between the two plots. In these 
models r^o = 0.35. 
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Fig. 3. — Poloidal energies at the core-mantle 
boundary for models with r^o = 0.35 (a), r^o = 
0.50 (b), and r^o = 0.70 (c) separated into non- 
axisymmetric components (upper) and axisym- 
metric components (lower). All points are a time 
average over at least one magnetic diffusion time. 
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Fig. 4. — Total poloidal energy at the top of the 
conducting layer {vm) as a function of mantle con- 
ductivity. All points are a time average over at 
least one magnetic diffusion time and have been 
normalized to the insulating mantle case. 



